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ABSTRACT 

Context. Circumstellar discs are expected to be the nursery of planets. Grain growth within such discs is the first step 
in the planet formation process. The Bok globule CB 26 harbours such a young disc. 

Aims. We present a detailed model of the edge-on circumstellar disc and its envelope in the Bok globule CB 26. 
Methods. The model is based on HST near-infrared maps in the I, J, H, and K bands, OVRO and SMA radio maps at 
1.1mm, 1.3 mm and 2.7 mm, and the spectral energy distribution (SED) from 0.9 /xm to 3 mm. New photometric and 
spectroscopic data from the Spitzer Space Telescope and the Caltech Submilimeter Observatory have been obtained 
and are part of our analysis. Using the self-consistent radiative transfer code MC3D, the model we construct is able to 
discriminate parameter sets and dust properties of both its parts, namely envelope and disc. 

Results. We find that the disc has an inner hole with a radius of 45 =b 5 AU. Based on a dust model including silicate 
and graphite the maximum grain size needed to reproduce the spectral millimetre index is 2.5 //m. Features seen in the 
near- infrared images, dominated by scattered light, can be described as a result of a rotating envelope. 
Conclusions. Successful employment of ISM dust in both the disc and envelope hint that grain growth may not yet 
play a significant role for the appearance of this system. A larger inner hole gives rise to the assumption that CB 26 is 
a circumbinary disc. 

Key words, circumstellar matter - planetary systems: proto-planetary discs - radiative transfer - stars: formation, 
individual: CB 26 



1. Introduction 

CB 26 is a small cometary-shaped Bok globule lo- 
cated about 10° north of the Taurus/ Auriga dark 
cloud at a distance of 140 pc (jLaunhardt fc SargentI 
|2QQ1[ ). An IRAS point source (IRAS 04559+5200) 
at its southwest rim su^ ^ests an embedded Class I 
young; stellar object (YSO) IStecklum et "aTl l2QQ4l ) source. 
iLaunhardt fc Henningj (|l997l ) found an unresolved 1.3mm 
continuum source associated with the IRAS source. 

Int erferometric observations by ILaunhardt fc SargentI 
(|2QQll ) showed that the major fraction of thermal dust emis- 
sion at millimetre wavelengths has its origin in a young 
circumstellar disc with a diameter of about 400 AU and a 
mass of about 0.1 M©. This disc is seen almost edge-on and 
the central star is not visible directly. However, the spec- 
tral energy distribution suggests a Class I YSO with L > 



0.5 L0 (jStecklum et al.ll2QQi ). From the ^^CO line e mission 
and t he Keplerian rotation curve, ILaunhardt fc SargentI 
(2001) derive a central stellar mass of = 0.5 ±0.1 M©. 
Furthermore, IL aunhardt et al.l (|2QQ8f ) detected a jet-like 
molecular outflow emanating perpendicular to the plane 
of the disc. This outflow seems to be co-rotating with the 
disc. 

We present new observations and a model for this source 
that accounts for spatially resolved data sets over more than 
3 orders of magnitude as well as for the unresolved SED of 
the object. 

Our modelling is based on spatially resolved maps of 
CB 26 in the millimetre regime from the Sub-millimetre 
Array (SMA) and the Owens Valley Radio Observatory 
(OVRO), high resolution images in the I,J,H7 K bands 
obtained with the the Near Infrared Camera and Multi- 
Object Spectrometer (NICMOS) at the Hubble Space 
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Telescope (HST) and the Advanced Camera for Surveys 
(ACS), which is as well an instrument of the HST. The 
photometric data for the SED upon our model is based 
are provided by the Multiband Imaging Photometer for 
Spitzer (MIPS) and the Infrared Array Camera (IRAC) 
aboard the Spitzer Space Telescope (SST) and millimetre 
photometry. We also obtained a spectrum with the Infra- 
Red Spectrograph (IRS) aboard the SST. 

In this paper, we use all available NIR to mm contin- 
uum data on CB26 to develop a self-consistent model of 
the source. This model consists of two parts: First, an opti- 
cally thick dust disc which accounts for the dark lane seen 
in images obtained with the HST of the object and the 
significantly elongated intensity profiles in the mm range. 
The second component is an optically thin envelope that 
reproduces the scattered light nebulosity. 

While the millimetre observations are sensitive only 
to radiation being emitted from dust in the dense region 
within the disc, the near-infrared images are dominated by 
scattered stellar light from dust in the circumstellar enve- 
lope and the disc's upper optically thin layers, often referred 
to as the "disc surface" . These observations trace different 
physical processes in different regions of the circumstellar 
environment, but they are both strongly related to the dust 
properties in the system. 

Thus, we are in the position not only to model obser- 
vations on the common basis of one set of parameters, but 
we are also able to investigate whether the dust proper- 
ties are different in the disc and the envelope. This has 
been suggested by investigations of the dust evolution in 
circumstellar discs where dust grain growth alters the dust 
grain properties in the circumstellar disc quite consider- 
ably whilst it is of less importance in the low-density en- 
velope. Evidence for grain growth has been fo und, for in- 
stan ce, for the c ircumstellar discs IM Lupi dP inte et al? 
[20081 GG Tau (iDuchene et al.1 l2004l : iPinte et al . 2007), 
HH 30 (IWatson StapelfeldtJ 120041). I RAS 04302+2247 
(IWolf et al.1 I2003D . and VV Serpens (lAlonso-Albi et al.1 

In order to compare our model with the available obser- 
vation s, we use the self-consistent r adiative transfer code 
MC3D (|Wolf et al.l[T999l : IWol j [2003br ) in a parameter space 
study on the free parameters of the model. We aim at find- 
ing the best-fit model, which we define to be the model that 
reproduces certain predefined features among the observa- 
tional data (such as width of the dark lane; see below for 
further details) best. 



2. Observations and data reduction 

In the case of the disc in the Bok Globule CB 26, we are in 
the fortunate situation of having a large variety of observa- 
tional data at hand. This data includes not only the spectral 
energy distribution (SED) from 0.9 /im to 2.7 mm but also 
resolved maps in the near-infrared and in the millimetre 
regime. We will now briefly discuss those observations and 
the data reduction in this section. 



2.1. HST Imaging 

The NICMOS and ACS data were taken by the GEODEH 
team. An overview of the complete program, its objects. 



^ Group for Edge-On Disc Exploration 



and its objectives can be found in Padgett et al. (2009, in 
preparation). 

2.1.1. ACS 

CB 26 was observed with the Advanced Camera for Surveys 
Wide Field Channel on 2005 September 09. Two 1250sec 
exposures were made in the F814W filter (A = 0.80 /im, 
A A = 0.15 /im), corresponding to Johnson I band. With 
a pixel scale of 0.05'^ the ACS provides a field of view of 
20f X 100'^ The images were reduced, combined to reject 
cosmic rays, and corrected for geometric distortion by the 
STScI pipeline. Residual hot pixels and cosmic rays were 
manually removed by replacing the affected pixels by local 
median values. 



2.1.2. NICMOS 

CB 26 was observed using NICMOS on 2005 September 15 
with the FllOW (A = 1.12 /im, AA = 0.16 /im ), F160W 
(A = 1.60 /im, AA = 0.12 /im ) and F205W (A = 2.06 /im, 
AA = 0.18 /im) filters on the NIC2 array. With a pixel 
scale of 0.075'^ the NIC2 array provides a field of view of 
19.2'' on a side. All data were taken in a MULTIACCUM 
step=128 sequence, with total exposure times of 512secs for 
the FllOW data and 384secs for the F160W and F205W 
data. The observations were dithered back and forth by 
0.75'' for bad pixel removal. 

All three of the NICMOS data sets were taken less 
than 1900 sees after emergence from the South Atlantic 
Anomaly (SAA) and therefore suffer from significant cos- 
mic ray persistence. The data were re-reduced using the 
most recent calibration files available. The standard STScI 
CALNICA pipeline was used as a basis for the re-reduction, 
in addition to which, we employed both the biaseq and ped- 
suh IRAF routines to remove the pedestal effects visible in 
each of the quadrants. The biaseq task removes the time- 
dependent variations in the bias level for each quadrant, 
and the pedsub task then removes the fixed bias offset. After 
the data have had the pedestal effect removed, we used the 
IDL S A ACLEAN piogidim to model and iteratively remove 
the persistent cosmic rays in the image. After emerging 
from the SAA, NICMOS powers back up and takes 2 dark 
frames. These dark frames provide the model for the cos- 
mic ray persistence pattern; SAACLEAN tokes this model 
and iteratively removes it from the data until the noise in 
the background reaches a global mi nimum. Further details 
on this procedure can be found in iBergeron fc Dickinson! 
(2003) (ISR 2003-010). Following the SAACLEAN pioce- 
dure, the data are run through a second round of pedestal 
subtraction, have any remaining bad pixels removed, and 
are then run through the CALNICB pipeline procedure to 
generate the final mosaics. 



2.2. Spitzer photometry and spectroscopy 

Mid- and far-infrared photometry for CB 26 were retrieved 
from the Spitzer Archive. IRAC observations were carried 
out under GTO program 94 (PI: Charles Lawrence). The 
data were taken 2004 February 11, AOR key 4916224, us- 
ing the high dynamic range mode and 5 dithered expo- 
sures of 30sec each. Aperture photometry was measured 
from the Spitzer pipeline "post-BCD" mosaics (version 
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11.0.2), using a standard 10 pixel radius aperture for the 
source and background annulus with 12-20 pixel radius. 
MIPS observations were made under GTO program 53 
(PI: George Rieke) on 2005 March 08, AOR key 12020480. 
A medium scan map covering 0.2° x 0.5° was performed, 
providing multiple dithered 4sec exposures and total inte- 
gration times of 168, 84, and 16 sec at 24, 70, and 160 
/im, respectively. Spitzer post-BCD pipeline mosaics (ver- 
sion 14.4.0) were used for aperture photometry. Photometry 
was measured in apertures whose radii and background an- 
nuh were 7^77'' - 13'' (24/im), 35'739'' - 65'' (70/im), and 
48"/64" — 128" (160 /im) with aperture correction factors 
as given by Engelbracht et al. ( 2007 ): Gordon et all (|2Q07) ; 
IStansberrv et all (|2007f ) 

We carried out our own Spitzer low-resolution spec- 
troscopy of CB 26 on 2006 October 18 under program 
30765, AOR key 18964992. The ramp duration and number 
of observing cycles used were 6x14 sec, 4x14 sec, 8x30 sec, 
and 8x30 sec for the SL2, SLl, LL2, and LLl modules of 
the IRS (respectively). The spectra were processed begin- 
ning with the intermediate droopres products from pipeline 
version 15.3.0. The 2D images were co- added and the nod 
pairs subtracted to remove stray light and other additive 
artifacts, as well as rogue pixels. ID spectra of each nodding 
beam were then extracted using a fixed aperture of 5 pix- 
els. Following Bouwman et al. (2008) the spectra were flux 
calibrated using a spectral response function determined by 
comparing standard star observations with template spec- 
tra. These response f unctions are also used with the "C2D" 
Legacy data set (e.g.'Kessler-Silacci et al. "2006). It was not 
necessary to apply any scaling to match the short-low and 
long-low modules, indicating that the source is unresolved 
at the IRS wavelengths. 

2.3. Millimetre and Sub-millimetre measurements 

2.3.1. SMA 

Observations at 270 GHz (1 .1 mm) with the Sub-Millimetre 
Array (SMA, Ho et al. 2004) were made in December 2006, 
in two configurations providing baselines in the range 12- 
62 kA. Typical system temperatures were 350-500 K. The 
quasar 3C279 was used for bandpass calibration, and the 
quasars B0355+508 and 3C111 for gain calibration. Uranus 
was used for absolute flux calibration, which is accurate 
to 20-30%. The data were calibrated using th e IDL MIR 
package (Qi 2005) and imaged using MIRIRAD ISault et all 
(|199 5'). The cleaned and restored 1.1mm continuum map 
was constructed with robust liv-weighting using line-free 
channels in both sidebands. Here we only use the contin- 
uum map. The observations together with the molecular 
line data are described in more detail in a forthcoming pa- 
per (Launhardt et al. in prep.). 

2.3.2. OVRO 

CB 26 was also observed with the Owens Valley Radio 
Observatory (OVRO) between January 2000 and December 
2001. Four configurations of the six 10.4 m antennas pro- 
vided baselines in the range 6 -180 k A at 2.7 mm (110 GHz) 
and 12-400kA at 1.3 mm (232 GHz). Average SSB sys- 
tem temperatures of the SIS receivers were 300- 400 K at 
110 GHz and 300- 600 K at 232 GHz. The raw data were 
calibrated and edited using the MMA software package 



(IScoville et al.l [19931 ). Mapping and data analysis used the 
MIRIAD toolbox. Observing parameters are described in 
detail in Launhardt & Sargent (2001). The data presented 
here include additional observations conducted in 2001. All 
maps are generated with robust i^'U- weighting, cleared, and 
restored with a clean beam. Effective synthesised beam 
sizes of all interferometric millimetre continuum maps used 
here are summarised in Table 13.11 



2.3.3. CSC 

Submillimeter observations of CB26 at 350 jam were ob- 
tained with the Submillimeter High Angular Resolution 
Camera II (SHARC-II) at the Caltech Submillimeter 
Observatory (CSO) on 2007 October 21. SHARC-II is a 
12 X 32 element bolometer a rray giving a 2.59' x 0.97' field 
of view (jDowell et al.ll20Q3[ ). The beam -size at 350 /im is 
8.5". 

We used the Lissajous observing mode to map a re- 
gion approximately 1' x 0.5', centred at the position of 
the source. We obtained two scans, each 10 minutes long, 
for a total integration time of 20 minutes in good weather 
(t225GHz ^ 0.06). During both scans the Dish Surface 
Optimisation System (DSOS j3 was used to correct the dish 
surface for gravitational deformations as the dish moves in 
elevation. 

The raw scans were reduced with version 1.61 of the 
Comprehensive Reduction Utility for SHARC-II (CRUSH), 
a publicly available^ Java-based software package. CRUSH 
iteratively solves a series of models that attempt to re- 
produce the observations, takin g into account both instru- 
mental and atmospheric effects (|Kovac s 2006; Kovacs et al.l 
[2OO6I : IBeelen et al.ll20 06). Pointing corrections to each scan 
were applied in reduction based on a publicly availably 
model fit to all available pointing data. Pixels at the edges 
of the map with a total integration time less than 25% 
of the maximum were removed to compensate for the in- 
creased noise in these pixels. We then used Starlink's stats 
package to assess the rms noise of the map, calculated using 
all pixels in the off-source regions. The final map has a la 
rms noise of 70 mJy beam~^. 

Photometry was measured in a 20" aperture centred at 
the peak position of the source. Calib ration was pe r forme d 
according to th e method used by Shirley et al.' (2000"); 
IWu et all (|2007l ). This method is based on the requirement 
that a point source should have the same flux density in all 
apertures with diameters greater than the beam FWHM 
(8.5" for these observations). To briefly summarise, a flux 
conversion factor (FCF) is calculated for a 20" aperture 
by dividing the total flux density of a calibration source in 
Jy by the calculated flux density in the native instrument 
units of juV in a 20" aperture. Flux densities of science tar- 
gets are then derived by multiplying the 20" aperture flux 
density (in the instrument units) of the source by the FCF. 
We measure a 350 /im flux density for CB 26 of 2.8 ± 0.6 

Jy- 



^ See www . CSO . caltech . edu/ dsos/DSOS_MLeong . html 
^ See www . subinin . caltech . edu/~sharc/ crush/ index . htm 
^ See www . submm . caltech . edu/~sharc/analysis/pmodel/ 
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Fig. 1. Inverse HST images obtained with ACS and 
NICMOS in the I, J, H and K band. The colour scale is 

2 

^ All four images have been rotated by —30° in order 
to align the major axis of the dark lane with the horizontal 
axis. For the image scale a distance to the object of 140 pc 
is assumed (100 AU = 0.7'^). The image orientation is PA 
330° up, and PA 60° to the left. 



3. Results from Observations - Basis for modelling 

The focus of this section are the results by the previously 
described observations. The presented data set forms the 
basis of our modelling. 



3.1. Images 

The maps from the HST in I, J, H, and K bands are shown 
in Fig. [H All all four images have been rotated in order to 
align the major axis of the dark lane with the horizontal 
axis. The emission seen on these maps is pure scattered light 
from the central star. The dark shadowy line that intersects 
the bipolar structure is present in all images. Especially the 
dependency of its width on the wavelength is clearly visi- 
ble. Since one expects the circumstellar disc to be optically 
thick at these wavelengths, the dark lane is interpreted as 
the disc's shadow in the encompassing envelope structure. 
The bipolar nebula also shows a complex morphology far 
above and especially below the disc. For further discussion 
we refer the interested reader to Padgett et al. (2009, in 
preparation). 

The interferometric millimetre continuum maps at 
1.1mm, 1.3 mm, and 2.7 mm, together with the dirty beam 
maps (see also Table [XT]) . are shown in Fig.O For the ease 
of modelling, these images have been rotated 30° in order 
to align the major axis of the elongated structure with the 
horizontal axis. While the source is unresolved vertically 
in all images, it is well-resolved along its horizontal axis, 
especially in the highest-resolution map at 1.3 mm. The 
1.3 mm map also recovers some extended emission from 
the envelope that might be related to a disc wind (see 






Fig. 2. Inverse, reconstructed images from millimetre in- 
terferometric observations {left column) and corresponding 
dirty beam maps {right column)^ linear colour scale. All im- 
ages have been rotated by —30° in order to align the major 
axis of the brightness distribution with the horizontal axis. 
For the image scale a distance to the object of 140 pc is 
assumed (100 AU = 0.7''). The contour lines are drawn at 
(from inside out) 90%, 50%, and 25% of the image maxi- 
mum flux value. In the dirty beam images, the 50% contour, 
marking the FWHM of the Gaussian clean beam, is marked 
bold. 



iLaunhardt fc SargentI [2QQll ) . For all images shown, radio 
and NIR, North is the same direction. 

Fig.[3]shows an overlay of a NIR colour-composite image 
and the 1.3 mm dust continuum emission. The spatial co- 
location of the millimetre dust emission and the dark lane 
in the scattered light images confirms the hypothesis of an 
edge-on optically thick disc as explanation for the observed 
features. However, the HST pointing is only good to about 
1 arcsec. Due to the high extinction in the cloud core and 
the small NICMOS field of view, there are no reference stars 
in the image that could be used to correct for this pointing 
uncertainty. 

3.2. Spectral energy distribution 

The results of photometric measurements are presented in 
Table 13.21 and FigHl The spectrum we obtained from the 



Table 1. Overview of the beam sizes in the minimetre-maps 



Instrument 


A [mm] 


PSF (FWHM) ["] 


Orientation 


SMA 


1.1 


1.00 X 0.84 


-59.3° 


OVRO 


1.3 


0.61 X 0.36 


-53.8° 


OVRO 


2.7 


1.21 X 0.87 


-58.3° 




Fig. 3. Overlay of the OVRO l.Srara continuuni raap and 
the HST NICMOS images. The NICMOS colour image is 
a 3-colour composite of F205W, F160 W, and FllOW (in 
the RGB colour planes, respectively), shown in log stretch. 
The contour levels are linear, with the lowest contour at 
about the 2.5 a-level, the others at the 6, 10, 14, 22, and 
25 cr-level, respectively. 



IRS very well complements the photometric points as seen 
in Fig. [11 

As is evident in this Figure, the 350 fim SHARC-II 
flux density is lower than expected based on comparison 
to the complete SED. This discrepancy can be explained 
by the fact that the Lissajous obser ving mode is insensitive 
to extended emission, as noted by Wu et al.l (j2QQ7lV This 
issue will be explored in more detail by M. Dunham et al. 
(2009, in preparation), but preliminary results suggest that 
flux densities measured from this particular observing mode 
may underestimate the true flux density by up to a factor 
of 2, even for relatively compact objects. 



4. Model & Modelling 

In this section we provide the reader with an introduction 
to the concepts and techniques we use. 



4.1. The Model 

First, we discuss the various components of our model. That 
is, the disc and envelope dust density distribution. The em- 
ployment of both these parts in the model is readily sug- 
gested by the data. 



in the Bok globule CB 26 
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Fig. 4. Spectral energy distribution. Data with error-bars 
from Table 13. 2[ The IRS spectrum is the solid line. The 
dashed-dotted line corresponds to the best-flt model. The 
dashed line is the best-flt model with the dust screen. 



4.1.1. The Disc 

The main part of our model for CB 26 is a parametrised 
disc. The disc is seen as the radially extended luminous 
structure in the millimetre maps. The dark lane in the near- 
infrared maps is another indication for a disc as it can be 
understood as the disc's shadow on the surrounding enve- 
lope. 

We employ a parametric approach to such a disc which 
can be written as 



Pdisc{r) = po ( — j exp 





' Z' 







(1) 



where z is the usual cylindrical coordinate with z = corre- 
sponding to the disc mid plane and rcvi is the radial distance 
from this z-axis (see e.g. lWolf et alll2003l : IStapelfeldt et aD 
'1998!). In our model, the parameter po is determined by the 
mass of the entire disc. is the stellar radius and /i, the 
vertical scale height, is a function of rcyi 



^(^cyl) = ^0 



rcyiV 



(2) 



Here the quantities a, /3, and ho in Equation ([2]) are ge- 
ometrical parameters. These parameters allow us to ad- 
just the disc structure and shape in order to flt the data. 
This modelling strategy has already been successfully ap- 
plied to various other edge-on se en discs, such a s the 
Butt erfly-Star IRAS 4302+2247 (IWolf et al.ll2003l). HK 
Tau (jStapelfel dt et al.|[T998i). IM Lu pi (jPinte et al.ll2008f ). 
and HV Tau (Stapelf eldt et al.ll2003[ ). 

Integrating Equation ^ along the z axis yields the sur- 
face density 



S(r) 



rcy\) 



(3) 



Comparison with Equation ([T]) yields the following relation 
between the exponent of the surface density power law and 
the geometrical parameters of the ansatz we use 



p = —/3 + a. 



(4) 
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A[/imJ 


Flux [mjy] 


Aperture [ J 


Instrument 


Keierence 


0.90 


0.062 ± 0.019 


24 


A TT A O r 

CAHA 3.5m 


/i \ 
(1) 


1.25 


2.2 it 0.2 


12 


A TT A O CT 

CAHA 3.5m 


/I \ 
(1) 


1.65 


8.2 ± 0.8 


12 


A TT A O r 

CAHA 3.5m 


/I \ 
(1) 


2.20 


17.1 =b 1.7 


12 


A TT A O C 

CAHA 3.5m 


/I \ 
(1) 


3.6 


18.3 ± 0.8 


12 


O 'U- TT~) A 1 

bpitzer IHACl 


4.5 


17.0 ± 0.8 


12 


O 'U- TT~) A /^O 

bpitzer 1HAC2 




5.8 


12.5 ± 0.7 


12 


O 'U- TT~) A /^O 

bpitzer 1HAC3 




8.0 


6.8 it 0.5 


12 


O '-i- TT) A /I 

bpitzer 1KAC4 




24 


160.6 ± 5.2 


30 


bpitzer MlPbl 




60 


A oor\ 1 orkr\ /i 

4880 ± 390.4 


75 


TT^ A O T)0/^ 

IKAb FbC 




70 


5555 =b 52 


120 


bpitzer JvllPb2 




100 


111 r»r\ 1 1 1 1 r\ 

11100 ± 1110 


125 


TT^ A O T)0/^ 

IKAb FbC 




160 


10731 ± 69.49 


24 


bpitzer MlPb3 




350 


2650 ± 850 


20 


CbO 






U 1 UU ZL lOUU 






l^J 


850 


600 =b 120 


54 


SCUBA 


(2) 


1110 


225 ± 45 


10 


SMA 


(2) 


1270 


240 ± 20 


54 


IRAM 30m 


(2) 


1300 


190 ± 30 


10 


OVRO 


(3) 


2700 


20 ±2 


10 


OVRO 


(3) 



Table 2. Photometric data points for CB 26. References: (1) IStecklum et af] (|2004l ): (2) Launhardt et al., in prep.; (3) 
iLaunhardt fc SargentI (|200lf ) 



The radial size of the disc rout is another parameter and 
is mainly determined by t he size of the e l ongat ed structure 
in the millimetre maps. As lHughes et all (|2008[ ) argues, the 
millimetre continuum does not trace the outermost region 
of the disc. However, in the case of CB 26, the extent of the 
dark lane is also consistent with the outer radius estimate 
we get. 

4.1.2. The Envelope 

In order to reproduce the pattern of scattered light we see 
in the observations in the I, J, H, and K bands, we add an 
envelope-like dust distribution to the model. 

The density within this envelope is understood to be 
orders of magnitude lower than that of the disc. The HST 
optical and NIR data show while the envelope has a high 
enough density to produce scattered light, it is orders of 
magnitude lower than the density of the circumstellar disc 
which is optically thick enough to obscure the central star 
completely. Further evidence for the low density in the en- 
velope are the millimetre maps where the envelope is cannot 
be seen. 

For the mo d el of the envelope structure we follow the 
ideas of lUlrichI (|l976l ). We thus implement a model for a 
rotating envelo pe resulting from in-fa l ling matter the sam e 
way as done bv lWhitnev et all (|2003f ): lEisner et all (|2005h : 

Here M is the dust in- fall rate, the stellar mass, the 
centrifugal radius and /i = cos 6. The initial in- fall path of 
dust particles is given by /io as r ^ 00. As this is the only 
occurrence of M and M^, the factor M(87r>/GM*r^)~^ in 
Equation [5] is merely a coupling constant that scales the 
mass of the envelope just as po does in case of Equation [H 
Hence, we are using M{S7r\^GM^r^)~^ = p as a fitting 
parameter. To avoid introducing a constant like p with no 
direct physical interpretation into the discussion later on, 
we shall in this work refer to and M separately. 



As a criterion to decide whether the disc dust distri- 
bution or the envelope dust distribution at a given point 
should be considered, we compare the two densities and 
choose the larger value: 

p(r) = I ^disclr-) : Pdiscir) > Penv{r) ^ .^x 

I Penv(''^) • Pdisc 

In this manner, we embed the disc into the envelope and 
guarantee a smooth transition from the disc to the enve- 
lope without the need to alter the density structure of the 
optically thick, millimetre-glowing part of the object. For 
the radius of the complete model space we take twice the 
outer radius of the disc. Since we do employ a maximum 
size for the outer disc radius rout, the remaining space from 
the disc's edge to the end of the model space is readily filled 
by the envelope. With a computational domain going out 
to 2rout7 we are able to model the scattered light from the 
envelope. 

4.1.3. The Dust 

Since gas is optically thirj3 in the wavelength regime we deal 
with, we limit ourselves to radiative transfer through the 
dust. For the mass relation of dust and gas we assume the 
standard value of = 100 which is in agreement with 

the findings of l OWe; et alJ in another disc sur- 

rounding a low- mass T Tauri star. Therefore, it is the dust 
whose density structure is described by Equations ^ and 
([5]) in the disc and the envelope, respectively. The dust grain 
properties in our model can be divided in three groups: The 
shape of the dust grains, their chemical composition, and 
their size distribution. 



^ This implies that we neglect line emission and absorption by 
the gas. We do not aim at modelling those as almost throughout 
the entire disc, the dust is by far the dominant coolant of the 
disc. Hence, thermal equilibrium obtained in radiative transfer 
calculations based on dust only will provide a reliable description 
of the disc's thermal structure. 
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Grain shape We assume the dust grains to be homogeneous 
spheres. Real dust grains, of course, are expected to feature 
a niuch more complex and fractal structure. As discussed 
by IVoshchinnikovl (|2QQ2[ ). chemical composition, size and 
shape of dust grains cannot be determined separately, but 
only as a combination. We therefore limit our model to the 
less complex but also less ambiguous approach of spherical, 
non-aligned and non-orientated dust grains. 

Grain chemistry For the chemical composition of the dust 
grains we employ a model that incorporates both silicate 
and graphite material. This grain model has alread y bee n 
used to model the "Butterfly star" by Wolf et all ()2QQ3[ ). 
For the optical data we use the complex refractive indices 
of "smoo thed astronomical silicate" an d graphite as pub- 
lished by 'Weingar tner fc Draind (|2QQlf ). Since the longest 
wavelength considered in our modelling is 2.7 mm, we ex- 
trapolate the refractive indices to that wavelength. This 
is readily done since for this wavelength regime both the 
real and the imaginary part of the refractive index show 
asymptotic behaviour. For graphite we adopt the common 
"I — I" approximation. That means, if Qext is the extinc- 
tion efficiency factor, then 

1 2 

Qext, graph = ^Qext(e||) + -^Qext{^±)j (7) 

where ey and e± are the graphite dielectric tensor's com- 
ponents for the electric field parallel and orthogonal to 
th e crystallo^raphic ax i s, resp ectively. As has been shown 
bv iDraine fc Malhotral (|l993l ), this graphite model is suf- 
ficient for extinction curve modelling. Applying an abun- 
dance ratio from silicate to graphite of 1 x 10~^^cm^H~^ : 
1.69 X 10~^^cm^H~^, we get relative abundances of 62.5% 
for astronomical silicate and 37.5% graphite (|e|| and 

Grain sizes For the grain size distribution we assume a 
power law of the form 

n{a) da ~ a~^'^ da with amin < a < «max- (8) 

Here, a is the dust grain radius and n{a) the number of 
dust grains with a specific radius. For amin = 5nm and 
<^max = 250 nm this distribution becomes the commonly 
known MRN distribution of the interstellar medium by 
iMathis et al.l (|l977h . We choose those values as the starting 
point of the present study. 

To model the different grain sizes and chemical popula- 
tions, one has to consider an arbitrary number of separate 
dust grain sizes within a given interval [amin • <^max]- But 
the observables derived from radiative transfer considering 
each grain species separately, are close to the observables 
resulting from radiative transfer (RT) simulations based on 
weighte d mean dust grain parameters of the dust grain en- 
semble (IWoljl2QQ3aD . Thus, we use weighted mean values 
for the efficiencies factors, cross sections, albedo, and scat- 
tering matrix elements. For each dust grain ensemble, 1000 
logarithmically equidistantly distributed grain sizes within 
the interval [amin • <^max] have been taken into account for 
each chemical component in the averaging process. 

There are arguments that the grain size can not be gov- 
erned by a power law as in Equation ([8|). Furthermore, in 
the complex environment of a circumstellar disc, we expect 
dust settling and grain growth to make the grain size dis- 
tribution quite dependent on the location within the disc. 



Unfortunately, a consideration of grain sizes that takes into 
account the effects of e.g. dust settling requires more than 
just one parameter as Equation ([8]) does. Given the data we 
have, we are not able to disentangle these parameters in our 
study. Therefore, we assume that the power law distribu- 
tion ([8]) to be valid in the whole disc and envelope structure 
and only use the maximum grain size as a parameter in our 
modelling efforts. 

It will turn out that even with these simplifying assump- 
tions we are able to model observational data of the system, 
although different amax in the disc and envelop e have been 
found - but for other objects - in the past (e.g. IWolf et al] 
I2003D . 

We assume an average grain mass density of Pgrain — 
2.5gcm~^. This density does not have any influence on the 
optical properties of the dust, as they are governed by the 
chemical composition and the particle size of the grains. 
The average grain mass density merely controls, together 
with the disc mass and particle size, the number of dust 
grains in the disc. 

4.1.4. Heating Sources 

There are two sources of energy for the disc that need 
our attention. The disc can be heated by stellar radiation 
and/or accretion of in-falling matter. 

Accretion Heating Our model involves a parameter with the 
dimensionality of an "accretion rate" : M in Equation ([5]) as 
part of the description of the envelope structure. However, 
this quantity is really an in-fall rate within the envelope 
and may not necessarily at the same time describe mass 
accretion onto the star itself. But it is in general the latter 
mass flow that, as e.g. in FU Orionis like objects, accounts 
for significant contributions to the system's luminosity. For 
a T Tauri like system as CB 26 we thus neglect accretion 
as a major source of energy. As our study shows, M is 
rather small, around ~ lO~^M0yr~^. This is small com- 
pared to FU Orionis objects but average for T Tauri Stars 
and strongly supports our ansatz. 

Besides matter in-fall in the envelope, accretion in the 
disc mi ght be an im portant source of energy. As shown 
by Wolf et al.l (j2003l ) , the accretion luminosity from within 
the disc is about two of magnitude smaller than the stellar 
contribution. As the model setup here is similar, we neglect 
accretion heating as a significant source of energy. 

Stellar heating The discussion above leaves the star as the 
only primary source of energy in our model. Its radiation 
heats the dust which then in turn itself re-emits at longer 
wavelengths. In this sense the disc in our model is passive. 
That is, we neglect accretion or turbulent processes within 
the disc as a possible other primary energy source. 

We do not observe the star directly. This has mainly 
two reasons: 

1. In the far infrared and at longer wavelengths, where the 
disc becomes less opaque with increasing wavelength, 
the contribution to the spectral energy distribution of 
the dust is orders of magnitude larger than of the star. 

2. In the optical, near, and mid-infrared bands, the disc 
becomes opaque and the star is shielded from our direct 
view. 
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Therefore, we have to assume the stellar parameters, that 
is temperature and luminosity, since we are not able to 
derive them directly from observations. Observations only 
hint at a luminosity being L > 0.5 Lq. As a starting 
point we choose an_^verage" T Tauri star as described by 
iGullbring et al.1 (|l998l ). This star has a radius of r = 2 
and a luminosity of = O.92L0. Assuming the star to 
be a black body radiator, this yields an effective surface 
temperature of Teff = 4000 K . 

Both of these parameters have been kept fixed in our 
parameter study to avoid degeneracies between parameters 
of the model. Except for the total flux, this choice has no 
impact on the near-infrared images. 

As lNattal (| 19931 ) showed, under certain conditions stellar 
light scattered back to the disc can have significant implica- 
tions for the thermal structure of the disc. Here, the outer 
envelope regions 400 AU < renv ^ 1000 AU is shown to be 
of importance. An important assumption in the argumen- 
tation is the scattering phase function to be independent of 
the scattering angle. However, in our modelling framework 
the scattering phase function is highly asymmetrical and 
favours forward scattering by orders of magnitude. Thus, 
the amount of radiation scattered back to the disc from the 
envelope outside our model space, i. e. at distance larger 
then 400 AU can be neglected. 

4.2. Means of modelling 

In this section we discuss how we proceed with the aforesaid 
model. We discuss the free parameters of the model, their 
range, and the sampling of the resulting parameter space. 
Finally, we review the constraints imposed upon our model 
by the various observation and give a criterion for the best- 
fit model. 

4.2.1. Radiative Transfer 

For our continuum radiative transfer simula t ions we mad e 
use of the program MC3D (Wol f et al.lll999l : IWolfl l2003bl ). 
It is based on the Monte-Carlo method and solves the 
continuum radiative transfer problem self-consistently. 
It estimates the dust temperature distribution taking 
into account any heating sources, in our case the cen- 
tral star's radiation. It makes use of the temperature 
correc tion technique as described by iBjorkman &: WoodI 
(|200l[) . the absorption concept as introduced by iLucvl 
(1 19 99*) and the enfor ced sc attering scheme as proposed by 
ICas hwell & Everett (|l959l ). The optical properties of the 
dust grains (scattering, extinction and absorption cross sec- 
tions, scattering phase function) and their interaction with 
the radiation field is calculated using Mie theory. Multiple 
and anisotropic scattering is considered. The phase function 
is highly asymmetrical (e. g. at the peak of stellar emis- 
sion at A = 0.7 /im one has (cos(^scatter)) = 0.86), strongly 
favouring forward-scattering. 

In order to derive a spatially resolved dust tempera- 
ture distribution, the model space has to be subdivided 
into volume elements inside which a constant temperature 
is assumed. Both the symmetry of the density distribution 
and the density gradient distribution have to be taken into 
account. For the present study, we use a spherical model 
space, centred on the illuminating star and an equidistant 
subdivision of the model in the ^-direction, whilst a loga- 



rithmic radial scale is chosen in order to resolve the tem- 
perature gradient at the very dense inner region of the disc. 
The required spatial resolution at the disc inner radius rim 
of our model ranges from 10~^ AU up to 10~^ AU and every 
grid cell outwards is 1% larger than its next inner neigh- 
bour. 

The radiative transfer is simulated at 101 wavelengths. 
The first 100 wavelengths are logarithmically distributed 
in the wavelength range [Amin, Amax-i] = [50 nm, 2.0 mm]. 
The largest wavelength used is Amax = 2.7 mm. 

With MC3D we compute observables from the model. 
These observables are then compared to the observed data 
in the quest for the best-fit model. Namely, the quantities 
we derive with MC3D from the model are 

1. Images in the NIR, that is in the I, J, H, and K Band, 

2. Images in the millimetre regime at 1.1 mm, 1.3 mm and 
2.7mm, 

3. 101 points for the SED accordingly to the above wave- 
length distribution. 

4.2.2. Constraints from observations 

Facing the broad variety of available observational data, one 
has to point out what the main features are that we want to 
reproduce with our model. This also determines the criteria 
for the best-fit model. The case is simple for the spectral 
energy distribution. There we aim at reproducing the com- 
plete spectrum over three orders of magnitudes from the 
optical bands down to the millimetre regime. We can di- 
vide the maps of the disc in two major groups: the maps in 
the millimetre regime, and the maps in the near-infrared. 
Both groups trace different physical processes and different 
spatial regions of the object. 

Millimetre maps For resolved images, the issue is not as 
simple as for the SED. Our model is rotational symmetric 
and thus does not provide for any related asymmetry as 
seen in observations. 

The morphology of the millimetre maps has its origin 
in the dust that is heated by the star and re-emits light at 
those wavelengths. Although the images at 1.1 mm, 1.3 mm 
and 2.7 mm are rather simply structured they impose two 
major features that constrain our models. These are 

1. the peak flux and 

2. the spatial brightness distribution. 

Since in all three maps the beam size is larger or comparable 
to the vertical extent of the disc, we can not constrain the 
flux distribution on the z-axis, perpendicular to the disc 
mid plane. Any feature there is smoothed out by the beam. 
Therefore, we focus on reproducing the flux distribution 
along the midplane of the disc. All images are fitted in the 
image plane. 

Maps in the near-infrared The four images in the near- 
infrared show more structures and details than the millime- 
tre maps. Besides the disc appearing as dark lane in the 
near-infrared also a complex, wavelength-dependent mor- 
phology of the surrounding envelope is seen. 

Considering that the circumstellar disc CB 26 is located 
at the edge of a Bok globule, one realises that the environ- 
ment of the disc can account for the majority of the sub- 
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structure seen on these maps, yet the Bok Globule itself is 
not part of our model. Hence, we have to restrict ourselves 
to the following two points that we want to reproduce with 
our model: 

1. the dependence of the width of the dark dust lane on 
wavelength and 

2. the relative peak height of the brightness distribution 
above and below the dust lane. 

We restrict our modelling efforts to a simple envelope struc- 
ture and the above two points since it is hard to distinguish 
whether the appearance of the object in the observations 
is due to envelope or environmental structure. We also put 
no emphasis on the vertical width of the upper and lower 
lobe nor the exact morphology. 

4.2.3. Quality of the fit 

For each comparison between model and observation on the 
aforesaid points, we get an individual xf • The total of 
one model is then just the sum over all the individual x^s: 

xLai = I xIed + X- + ^ X • (9) 

\ mm— maps NIR— maps / 

Here, n = 8 as we have one x^ from the SED, three from the 
millimetre maps and four from the scattered light images. 

Based on Xtotai from Equation ([9]), we give our 

modelling errors as the range where we can alter the param- 
eter values without changing Xtotai ^oi"^ than 10%. This 
value is rather arbitrary as there is no mathematical rea- 
soning behind it. Yet, it has proved within our study to re- 
flect quite well the adjustability of the model. Allowing for 
a larger variation of Xtotai than 10% gives generally worse 
results. 

4.2.4. Parameter space study 

Based on the model outlaid in the previous sections, we 
are left with ten adjustable parameters to reproduce the 
characteristics as described in the previous section. If not 
stated otherwise, we choose the range of a parameter for our 
stu dy based on mod elling of other yet similar objects (see 
e.g. IWolf et al.ll2QQ3 ). Then we first sample that range of 
an individual parameter in four coarse steps, select the two 
best values and go with the same procedure to the next pa- 
rameter. Secondly, we take the results as an indicator how 
to refine the stepping in which smaller range. This process 
is being iterated to reach our final results. As the starting 
values in our parameter space study we chose the values 
obtain for circumsteUar disc IRAS 04302+2247 which has 
at first glance a similar appearance as the disc in CB 26. 
In detail the parameters we have are 

1. The exponents a and /3 which describe the radial density 
profile of the d i sc (se e Equation ([T]) for details). From 
[D'Alessi o et al.l (| 19991 ), we choose for the flaring param- 
eter /3 = 1.25 and get then a corresponding value for 
a = 2.25 from the relation a = 3( 0—^), which is a result 
of acc retion disc physics (see e.g. IShakura fc SvunvaevI 
(|l97 3hV Those values are taken as a starting point and 
we then look for agreement between observations and 
modelling at and beyond those values in steps of 0.1 
and 0.2, respectively. 




Fig. 5. Sketch showing all components of our model and 
their spatial arrangement. All sizes are not to scale. 



2. The scale height ho of the disc at a given radius. In the 
following we will fix the scale height at a radial distance 
from the star of Tcyi = 100 AU and consider the cases 
ho = 5 AU, 10 AU, 15 AU, 20 AU, and 25 AU. 

3. From the density Equation ([5]) for the envelope we have 
the centrifugal radius rcf and 

4. the coupling constant p which is a function of the 
mass accretion rate M and the mass of the central 
star M*. For the centrifugal radius we probe values 
in the range 100 AU < rcf < 800 AU. As for p, we 
chose for the accretion rate values from the inter- 
val M e [10-^ Moyr-^ 10-^0 Moyr-i] and calcu- 
late p under the consideration of the stellar mass by 
iLaunhardt et al.l (2008). In this paper a dynamical mass 
of M* = 0.5 ± 0.1 Mo is derived. 

5. The inner and outer disc radius. The inner radius was 
initially set to 0.1 AU, which is approximately the dust 
sublimation radius. However, we do not get good agree- 
ment especially with the 1.3 mm map unless we choose 
values of ~ 45 AU. For the outer radius, we chose a 
value of 200 AU but also considered configurations with 
150 AU and 250 AU. 

6. The disc mass. We consider 7 different disc dust masses 
in the range from 1.0 x 10"^ up to 4.5 x 10"^ Mq. 

7. The disc's inclination is such that the system is seen 
almost edge-on. We restrict the range for 6 to be part 
of our parameter space to values between 60° and 90° 
with a stepping of 1° between 80° and 90°. 

8. The maximum grain size amax- Grain growth is a ma- 
jor issue for protoplanetary discs as it is the first step 
towards the formation of planets from the dust in the 
interstellar medium. We allow for a maximum grain size 
of 1 mm. 



For a summary of those ten adjustable parameters of our 
model and the range we covered in the study see Table [3l 
A sketch showing all components of our model is presented 
in Fig. [51 
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Dust mass [Mq] 


1.0 X 10"^ 
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±0.2 X 10-^ 


M[M0yr-i] 
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Table 3. Overview of parameter ranges and best-fit values. For the definition of the uncertainty see section 3.2.3. The 
first group of parameters contains purely geometric parameters, the second group physical parameters and the last group 
the inclination of the disc as seen by the observer as an observational parameter. 



5. Results 

The values of the parameters of our best-fit model can be 
found in Table [3l Our geometrical parameters a and P of 
the disc density structure yield with the Equation (|4]) a 
surface density power-law exponent of p = —0.8. 

In Fig.[6]the density and temperature distribution of the 
best-fit model is shown while in Fig. [71 radial profiles of the 
density and temperature distribution in the midplane and 
20 AU above are shown. One can see, that the high density 
in the midplane on the inner rim provides enough opacity 
such that material behind it in the midplane is not being 
heated up directly by the stellar radiation. The highest dust 
temperature at the inner disc rim amounts to ~ 90K and is 
reached ~ 20 AU above the midplane. However, due to the 
high density in the midplane, the stellar radiation does not 
penetrate deep into the midplane which results in the steep 
temperature gradient. In the less dense upper (and lower) 
layers of the disc, the stellar radiation also heats more dis- 
tant parts of the disc resulting in a less steep temperature 
gradient. 

An average dust temperature of Tdust = 16 K is obtained 
from the temperature distribution by weighting it with the 
mass distribution. This goes nicely with Fig. [6] if one bears 
in mind that the bulk of the dust is located in the midplane 
and well shielded against stellar radiation by inner parts 
of the disc. High temperatures are only reached in a very 
narrow region at the inner disc rim and in the very low 
density regions of the disc and thus contributing little to 
the mass averaged dust temperature. 

Since our dust grain model only uses refractive indices, 
an effective dust grain opacity can be calculated by re- 
arranging 

Mdu.t = -^—-^ (10) 

and assuming gas-to-dust ratio of 100. Here, hCi, is the wave- 
length depended mass absorption coefficient, Sj^ is the ob- 
served flux, D the distance of the object, and (Tdust) 
the Planck-function at a certain temperature. The calcu- 
lation yields for our model a dust opacity of /^i.3mm = 
0.26cm^g~^. Th is value is very close to the ISM dust 
opacity given by iDraine Led ((l984l . ll987f ). Compared to 
opacities for co agulated dust grains and ice- c oated grains 
lOssenkopf fcHennin g (1994); Beckwith eTall |l990), our 
model yields an effective dust grain opacity at the lower 
end of the range of commonly employed opacities. 



Fig.Hlshows the spectral energy distribution of the best- 
fit model in comparison with the spectral data. We achieve 
quite a good match to the observational data except for the 
optical wavelengths. Since the circumstellar disc is embed- 
ded in the Bok globule CB 26 we need to care about the 
dust outside our model space as well. Thus, we add a screen 
that mimics the effect of foreground extinction between the 
object and the observer. For this screen, we assume the ex- 
tinction properties of inte rstellar dust gr a ins. S uch a screen 
is descriloed in detail bv ICardelli et al.l (|l989l ). Using Ay 
as a parameter with a minimum value of 2, we find in our 
study that a screen with a visual extinction of Ay = 10 can 
easily account for the missing flux in the optical. The result 
is shown in Fig. [Has the dashed curve, whereas the dashed- 
dotted line corresponds to the best-fit model without the 
screen. 

It needs to be stressed, that this screen only applies ex- 
tinction law to all observable quantities. It is not subject to 
any radiative transfer or thermal re-emission. We estimated 
the possible contribution to re-emitted radiation of a such 
a screen with Ay = 10 composed of ISM grains at the same 
distance as CB 26 with a temperature of 16 K. We found the 
screen to be clearly optically thin (e.g. ri.3mm = 8 x 10~^) 
and has only enough mass to have about ~ 1% of the ob- 
served flux in the millimetre regime. 

Furthermore, the spectral energy distribution of the 
model shows that the contribution of the envelope is quite 
important for shorter wavelengths. In this regime, the main 
contribution to the spectral energy distribution comes from 
the envelope whilst in the radio regime the flux comes com- 
pletely from the dust in the disc that glows at those wave- 
lengths. Fig. [8] illustrates this. 

6. Discussion 

The following paragraphs are now dedicated for a more 
detailed discussion of some results of our model. 

6.1. Grain size and growth 

It needs to be pointed out, that we found a model capable 
of explaining all major elements of the observations without 
the need to increase the maximum grain size in our parame- 
ter study significantly. The maximum grain size of our best 
fit model is amax = 2.5 /im. While this is a factor ten larger 
than the smallest maximum grain size considered in our 
parameter space, the value found is only marginally larger 
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Fig. 7. Radial profiles of density {upper plot) and tem- 
perature {lower plot) in the midplane along r /me) 
and along z/r = 0.1 {dashed line). The latter profile is cho- 
sen in order to include the scale height /iq at r = 100 AU. 
The maximum density is normalised to 1 (corresponding to 
0.34 X 10-^gcm-^). 
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Fig. 6. Upper plot: Contours of dust density distribution 
on a plane perpendicular to the midplane normalised to 
the peak density of 0.34 x 10~^gcm~^. The contour levels 
are at 10"^, 10"^, 10"^ 5 x 10-^ and 2 x 10"^ Lower plot: 
Contours of the temperature distribution on the same plane 
as above. Contour levels are at 20K, 40K, 60K, 70K, and 
80K. The maximum temperature is 90K. 



then upper grain sizes given for the ISM in the literature. 
In fact, this is only true for the dust in the the disc compo- 
nent of our model. The maximum grain size in the envelope 
is the same as in found in the ISM, amax = 0.25 /im. If the 
maximum grain size in the envelope were bigger, then the 
short wavelength part of the SED could not be reproduced. 

It is certainly smaller by several orders of magnitude 
than the maximum gr ain size found in o ther disc models 
such as in the work of iPinte et al.l (|2008[ ). There, a max- 
imum grain size of a few millimetres has been found. In 
contrast, models in our parameter space featuring values of 
<^max ^ 1 mm fail to fit the SED in the millimetre regime 
as the slope of the model SED is not steep enough. Also, 



Fig. 8. Contributions to the spectral energy distribution 
from disc (dashed line) and envelope (solid line). The tran- 
sition from envelope to disc as the major source of radiation 
(re-emission and scattering) is at A = 217 /im. 

we did not succeed to reproduce the default value for the 
maximum grain size of amax = 250 nm of ISM. In particu- 
lar, the model would be off by a factor of ten for the SED 
data point at 1.3 mm. 

One needs to discuss why the slight change for amax 
from 250 nm to 2.5 /im allows for a fit in the millimetre 
part of the SED - at wavelengths three orders of magnitude 
larger than the largest grains in the model. Intuitively, one 
would expect the millimetre part of the SED to remain un- 
altered by a change of grain size at that level. However, this 
expectation is based on the assumption that the absorption 
efficiency of the grains Cabs in the millimetre regime is also 
insensitive to a change in the grain size at the micrometre 
level. Yet, this only holds true for only two of the three 
dust species in the model. For the astronomical silicate and 
the graphite component with an alignment of the crystals' 
optical axis perpendicular to the propagation direction of 
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the electromagnetic field Cabs has the same slope in the 
millimetre regime for grain size distributions with a maxi- 
mum grain size of 250 nm and 2.5 /im. But for the third dust 
species, namely the graphite component with the crystals' 
optical axis aligned with the electromagnetic field, this is 
different. Here, the slope of Cabs is significantly larger for 
<^max = 250 nm than for amax = 2.5 /im. This effect is large 
enough to dominate the sensitivity of the SED to changes in 
the maximum grain size even at the level discussed despite 
the fact that the dust species responsible for this behaviour 
has only a 12.5% share of the total dust. This is due to 
the fact that graphite is a far more effective absorber than 
silicate. 

A look on the millimetre spectral index of the data 
yields amm = 3.1 ±.27. The corresponding millimetre opac- 
ity slope is Pmm = 1.1 ± .27, if the millimetre emission is 
assumed to be optically thin. A /3mm = 1 and smaller is 
understood to indicate dust grain particles larger than in 
the interstellar medium to be present. A value of Prnm = 2 
is expected if only ISM grains were present in the disc. The 
latter is true only for grains whose absorption efficiency 
Cabs behaves like silicate. The value of Prnm we obtain from 
data and model of CB 26 close to what is expected for large 
grains despite having still only micrometre sized grains in 
the model is due to the unorthodox behaviour of the par- 
allel graphite component on the one hand side and on the 
other hand side to the non- vanishing optical depth in the 
millimetre regime (e.g. ri.smm ^ 0.6). 

iDraind (^2006) considered also the behaviour of the mil- 
limetre opacity slope Prnm for dust mixtures of graphite and 
silicate. There, no dependency of the opacity index as in our 
work was found. However, in this analysis not crystalline 
graphite was used but the optical properties of amorphous 
carbonaceous solids instead. 

We summarise that we do not need mm sized grains 
to model the circumstellar disc CB 26 but grains with a 
maximum grain size still close to what is found in the ISM. 
This is in contrast to the modelling of the Butterfly star 
(see IWolf et aL 20031) as well as for the circumstellar disc 
HH 30 fsee lWood et al.ll2002[) where in both cases the au- 
thors found it necessary, to have their largest grains at least 
four times larger than the largest grain of the interstellar 
matter. Of course, this result is based directly on the choice 
of the grain model we made. For another model, especially 
one without graphite, larger grains might be needed to fit 
the observed SED. Yet, due to the poorly constrained dust 
composition of circumstellar discs - in particular in the disc 
interior - this degeneracy between dust model and grain 
size must remain. As the dust model used in this work is 
also used in the context of other studies of circumstellar 
discs such as the Butterfly star (|Wolf et al.ll2003f ), it is a 
reasonable choice as it keeps the models of similar objects 
comparable as they are built on common assumptions. 



6.2. Inner hole 

The most unexpected result of our modelling is the inner 
disc radius for the model. In Fig. [9] the solid line shows the 
flux profile along the disc midplane at 1.3 mm as seen by 
OVRO whilst the thin solid line is the PSF of the observa- 
tion. At the centre of the disc, the profile shows a plateau 
in the brightness distribution. There are two explanations 
in stock for this dip. 



Horizontal Cut 



3mm 




Fig. 9. Horizontal cut through the spatial brightness dis- 
tribution at 1.3 mm. The thick solid line represents the 
OVRO observation, the dash-dotted lines gives the addi- 
tion/subtraction of one <j, the dotted line indicates the 2a- 
levels, the dashed line corresponds to our best-fit model, 
the thin solid is the dirty beam of the observation. 



First, the minimum could indicate that emitting dust in 
that region is present but not visible. If the optical depth in 
the midplane is sufficiently high, the flux contribution from 
the inner parts of the disc compared to the contribution 
of the optically thin parts on the disc's surface would be 
smaller. 

A high optical depth can easily be reached in regions of 
high dust densities. For a given total disc mass, a disc with 
small inner radius yields higher densities than the model 
with the large inner radius. Hence, the smaller the inner 
radius in our model, the more matter we find to be close to 
the star and thus reaching higher optical depths in the inner 
disc regions. However, even for the smallest inner radius of 
our parameter space, 0.1 AU we did not reach an optical 
depth that obscures enough flux from the disc centre. 

This behaviour would also be more obvious when com- 
pared with maps at shorter wavelengths since the optical 
depth increases with decreasing wavelength. But in the im- 
ages at 1.1 mm and 2.7 mm we do not observe a dip at the 
centre of the disc. Unfortunately, the available images do 
not help us to conclude whether the absence of the dip is 
really an indicator for a big void in the disc. This is be- 
cause the point spread function at those two i mages is far 
too large to resolve the feature (see Table \3A\i . I Wolf et al.l 
( 2008, ) reasoned this way in the case of IRAS 04302-1-2247 
where they found a similar dip in the brightness distribu- 
tion at A = 894 /im but not at 1.3 mm. 

The second possible cause for the spatial brightness dis- 
tribution in the 1.3 mm map is the actual lack of dust in 
the inner region of the disc. Whilst we started with an inner 
radius in the order of magnitude of a few tens of the stellar 
radius, it was not possible to match the plateau structure 
of the 1.3 mm map. On this account, we allow the inner disc 
radius to be as large as 50 AU. Unfortunately, within our 
model and parameter framework we are not at liberty to 
increase the total mass, and hence the density at the disc 
centre because the flux at 2.7 mm sets already an upper 
limit for the total disc dust mass. This is because at this 
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wavelength any model is optically thin and thus we see the 
total matter of the disc. 

Another way to think about the issue is to consider the 
optical thickness of the disc at 1.3 mm. We have run one 
simulation of the disc with exactly the same parameter val- 
ues as for our best-fit model except for the inner radius. 
For the comparison model we chose rin = 0.1 AU. For both 
runs we computed the optical depth along the line of sight 
from the observer through the disc. As a result, for the 
large inner radius we have ri.smm = 0.6 and for the small 
inner radius ri.smm = 1-9- In the first case we deal with 
an optically thin system. As for the second, this is not so 
clearly said. An optical depth of 1.9 means that the initial 
flux is reduced by a factor of e~^-^ = 0.15. A much larger 
value of r along the line of sight would be required to hide 
all emitting dust in the central region and produce the ob- 
served plateau structure. In our study, we were not able to 
fit the millimetre profile with a small inner radius. 

We are therefore forced to conclude, that in the mil- 
limetre regime we do see the entire disc. Hence, the plateau 
structure rises from the wide spatial separation of the inner 
rim from the star, whereas a disc with a small inner radius 
would have a central peak in the brightness distribution. 

Based on this line of arguments, our model provides 
predictive power for high resolution images at wavelengths 
longer than 1.3 mm. Since we conclude that the plateau 
structure in the brightness profile is due to lack of dust, it 
should also be visible at longer wavelengths. There, the dust 
becomes ever optical thinner and thus provides us with lit- 
eral insight inside the disc. Now, if the lack of flux in the disc 
centre would be due to an effect of optical depth, it should 
vanish with longer wavelengths and the brightness profile 
should have one central peak instead of a plateau-structure. 
Unfortunately, the image at 2.7 mm has a point spread func- 
tion far too large to allow us to disentangle between the two 
possibilities. Future observations of the disc CB 26 with a 
high spatial resolution will provide a perfect opportunity 
to confirm this prediction. Fig. [10] shows what we expect 
the disc to look like at different wavelengths according to 
our model with the inner void. With the Atacama Large 
Millimeter Array (ALMA) it will be quite easy to confirm 
our findings. 

It needs to be pointed out, that the spectral energy 
distribution from CB 26 itself does not hint on the large 
inner hole. Its general shape is similar to other edge-on 
seen disc's. This is because the inner regions of the disc, 
while they completely dominate the disc emission in the 
1-20 micron regime, is completely hidden from view in the 
edge-on configuration. Scattered starlight from the outer 
disc and/or the envelope accounts for the short wavelength 
bump of the SED, irrespective of the amount of emission 
from the disc itself. Only the quality of the OVRO map at 
1.3 mm provides us with the possibility (and the need) to 
conclude about the inner clearing. 

A possible explanation for a large void with 90 AU in 
diameter as we found it here might be, that the disc in the 
Bok globule CB 26 is actually a circumbinary disc. Binarity 
is also a possible ex planation for the rotating molecular 
outflow described bv iLaunhardt et aH ([2008^ ). However, a 
detailed dynamical study is not the scope of this paper. 

Of course, another idea might be that the dust in the 
inner region has already bean processed to planetesimals, or 
at least to bodies that are large enough to decouple from the 
disc an its dynamics and do not contribute to the mm-flux. 
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Fig. 10. Predicted appearance of the inner region of the 
disc. 



However, a first indicator for a low age of the system is the 
fact that it is still deeply embedded in its parental cloud. 
Another indication to the young age of the system is the 
absence of dust grains larger than found in the interstellar 
medium. Hence, we do not expect that we may be dealing 
with a so called "transitional disc"0 and a planet population 
in the centre that has already cleaned up those formerly 
dusty regions. 

6.3. Disc mass 

Within our model a rather high disc mass is needed to 
account for the observed flux in the millimetre regime. With 
a total of about 3 x 10~^ Mq and dust-to-gas ratio of 



100, 



(11) 



we end up with a total disc mass of ~ 0.3 M© which is close 
to the star's mass. Hence, we need to reconsider the mass 
we computed for the disc. The derived disc mass essentially 
builds upon the assumptions we had to made about the dust 
grain chemistry and shape. 

6.3.1. The dust grain structure, temperature and density 

The radiative transfer in our modelling is performed under 
the assumption of spherical dust grains to avoid equivo- 
calities. A possibility to have a lower estimate of the dust 
mass is to dismiss this assumption. This very much compli- 
cates the radiative transfer calculations since we now need 
to take into account all the possible fractal grain shapes 
as well as the spatial orientation of every grain. So far, we 
do not have the abilities to do such a simulation. But we 
can make a gedankenexperiment on parts of the results. 



^ Those discs are considered to be predecessors of evolved 
debris-type discs. 
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Of course, one expects strong changes in the scattering be- 
haviour of dust grains. But besides that, one also can think 
of a plenitude of dust grains with almost the same absorp- 
tion cross section as spherical dust y ains but with much 
less mass. IVoshchinnikov et al.l (|2QQ7[ ) discusses very fluffy 
particles with a porosity up to 90%. 

We furthermore want to point out, that the grain den- 
sity Pgrain IS uot ouc of the fitting parameters. The disc's 
mass is proportional to this density and the number of 
grains. We only can constrain within our study the disc 
structure, the grain size and their number, but not the 
density of one grain. For our investigation of the system 
we used Pgrain = 2.5gcm~^, but it might well be less than 
this value. In turn, this will alter our estimate for the disc 
mass by the same factor. 

6.3.2. The snow line 

For the estimation of the total disc mass, gas and dust, one 
makes use of the canonical relation (pTj) . This relation as- 
sumes, that there is some gas and dust of the disc inside 
the snow line and some outside. The snow line indicates the 
largest radius for which the temperature in the disc is high 
enough to keep water from freezing onto the dust grains. 
This line is usually set at a radius where the disc temper- 
ature drops to 170 K. As reported above, the maximum 
temperature we reach within our disc is about 90 K, and 
the average is only 16 K. This means, that our entire disc 
is outside the snow line. Thus we need to ad just the dust to 
gas ratio (pT]) to about 50 as published bv iKokubo fc Idal 
1002.). 

However, the very same dust model we made use of in 
this work was a lso build upon in other modelling projects 
(eg. IWolf et al.L 2003). In order to allow comparison of our 
model for CB 26 with those models, we also give considera- 
tions about the disc stability when we keep the dust model. 

6.3.3. Stability, Binarity &l Dynamics 

A criterion for a disc to become unstable was first shown 
by iToomrd (|196 4) . In order to allow self-gravity of the disc 
to take over, it requires to satisfy 



where Cg is the local sound speed, n the angular frequency 
of the disc, G the gravitational constant and E the sur- 
face density. For values of Q smaller then unity, the disc 
is assumed to be gravitationally unstable whereas for Q 
larger unity the disc is supposed to be stable against grav- 
itational collapse. As Fig. [11] clearly shows, in our model 
Q <1 throughout the entire disc by a factor of two to five. 

Yet, the Toomre criterion only can be consistently em- 
ployed for a system that is formed by a central star and 
a surrounding disc. As discussed above, the inner hole we 
find in CB 26 can be an indication for a binary. An exam- 
ple for such a system would be the young binary system 
GO Tau. Here, a dusty ring around the central stars has 
been observed by iGuilloteau et al.l (|l999l ) with a total ring 
mass of 0.13 M0, which is about a factor of two small than 
the total mass of CB 26. Yet, having two stars and the disc, 
the discussion of stability of such a system is quite delicate 
and not the topic of the present paper. 
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Fig. 11. Contour lines of the Toomre parameter Q (see 
Equation [T2|). The lines, from right to left, are at levels of 
Q = 0.1, 0.3, 0.5, 0.7. 



However, if one accepts to hold on to a single central 
star, there are still effects that we did not take care of in 
our model but might be important for the system's stabil- 
ity. Following the discussion in the paper of lGammid (|200lf ). 
discs violating the Toomre criterion are not stable but nev- 
ertheless might be in a steady gravoturbulent state. The 
paper investigates gravitationally unstable thin Keplerian 
discs and concludes, that the actual outcome of the insta- 
bility depends on the cooling time r^. 

However, the parametric ansatz we use does not account 
for any dynamical interaction within the dust. Therefore, it 
would be interesting to couple our radiative transfer code 
MC3D with hydrodynamic simulations of circumstellar discs 
in a future study. 

6.4. The width of the dark lane &i envelope structure 

Another aim of our modelling efforts was to mimic the chro- 
maticity of dust lane, which is narrowing with increasing 
wavelength. Fig. [T2l shows cuts from north to south through 
the centre of the NIR maps from modelling and observa- 
tions. Our model of the spherical envelope is quite success- 
ful in reproducing the overall flux at each wavelength. Also, 
wavelength-dependent width of the dark lane is correctly 
reproduced. 

Fig. [13] shows an overlay of the HST/NICMOS H-band 
image with the contours given by the best-fit model. The 
Fig. clearly shows, that the dark lane is well reproduced. 
Since our envelope model. Equation ([5]), is axial symmetric, 
we cannot expect the asymmetries of the lower lobe to be 
modelled as well. Thus, the model yields the seen discrep- 
ancy to the observation. 

A second effect of the aforesaid disregard is that the 
lower lobe flux, which in the image is concentrated on the 
left hand side, is in the model distributed among the com- 
plete lower lobe. Since both fluxes are equal in magnitude 
we end up with a smaller maximum in the model which 
also can be seen in Fig. [13] 
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Fig. 12. Vertical cuts at the horizontal centre through the 
near-infrared scattered light images. The solid line is from 
observational data, the dashed line from our best-fit model. 
For an illustration where those cuts are obtained in the re- 
spective maps, see Fig. [131 The plots are normalised to 1. In 
case of the TBand, this corresponds to 1.3 x 10~^ Jy/beam, 
5 X 10~^ Jy/beam in the J-Band, 2 x 10~^ Jy/beam in the H- 
Band, and 2.2 x 10"^ Jy/beam ('beam' refers to the FWHM 
area of the PSF). 
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Fig. 13. Inverse H-Band image from HST/NICMOS. The 
contour lines are from our best-fit model. From the outer 
to the inner lines are they at 27%, 41%, 68%, of the peak 
height. The second dark lane is not reproduced by our mod- 
elling as we did not aim at envelope asymmetries. The PSF 
for those images is rather small ~ 7AU FWHM. The ver- 
tical line illustrates how the cuts were taken for Fig. [121 

6.4.1. Alternative model for the envelope 

Besides the rotating, in-falling envelope model as described 
in Equation ([5|), we also tested a spherical symmetric den- 
sity distribution as a model for the envelope: 

Penv{r) = pCenv^^yi (13) 

This brings two parameters into the model, po,env and 7. 
With this parameters we are able to model the dependency 
of the dust lane wideness on wavelength as well as the over- 
all SED. Yet, this approach has two major drawbacks: 
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Fig. 14. J Band image from a simulation with a power-law 
envelope structure (left) and from rotating envelope (right). 
See also Fig. [H for the observation. 

Model @ 1 .3 rnm — Different inclinations 
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Fig. 15. Dependency of the brightness distribution along 
a central cut on the inclination. All the curves for = 
90°, 85°, 80° overlay. 

1. The spatial flux distribution as seen in the I, J, H, and K 
images cannot be reproduced. In fact, we obtain a much 
more concentrated flux distribution above and below the 
dark lane than what we see in Fig.[3l Fig. [T4l illustrates 
this. 

2. In the SED appears a strong silicate emission feature 
between 8 jam < A < 10 /im. We are not able to have the 
feature disappearing as required by observations except 
for the inclination approaching values 75° > 0. But this 
clearly contradicts the edge-on nature of the system. 

Therefore, we discarded Equation ([T3]) as a model for the 
envelope. 

6.5. Inclination of the disc 

Within our study, we considered two ways to get a hand on 
the disc's inclination. The first one is to infer it from the 
millimetre observations and the second one from the maps 
in the near-infrared. 

In Fig. [151 the spatial brightness distribution along a 
horizontal cut in the 1.3 mm map from the best-fit model 
at different inclinations is shown. For inclination values in 
the range of = [80°; 90°] the profile does not change con- 
siderably. For inclinations smaller than 80° the profile gets 
lower. This is because for an observer the warm inner rim of 
the disc looks at exactly 90° like a line and becomes more 
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and more a circle as one goes from edge-on to face-on ori- 
entation. As long as the complete rim fits into the size of 
the point spread function, one cannot distinguish between 
the different inclinations. But as soon as the rim gets less 
ellipsoidal and is not longer within the scope of one beam, 
the flux in the centre decreases. This is what can be seen in 
Fig. [151 Thus, from the comparison of models at different 
inclinations with the observations, we can infer an value for 
theta in the aforesaid range of ^ = [90°; 80°]. Naturally, this 
way of reasoning is only valid for optically thin systems. As 
outlined earlier, this is the case for CB 26. 

The second possibility to conclude about the disc's in- 
clination is to compare the relative peak heights from verti- 
cal cuts through the spatial brightness distributions in the 
near-infrared. Assuming a symmetric system one expects 
that both peaks exhibit the same height if the inclination is 
exactly edge-on. From an analysis of the actual peak heights 
(see Fig. [12]) in the HST/NICMOS images we deduce an in- 
clination of ^ = 85° ± 5° which is in good agreement with 
the numbers we obtained from the millimetre maps. 

6. 6. Comparison with similar objects 

In general our model has disc and envelope parameters that 
are comparable with those of the circumstellar environment 
of other young stellar objects. Yet, there are not many edge- 
on seen circumstellar discs in the sky for which the same 
richness of observational data is available for modelling as 
it is for CB 26. Hence, the number of comparable multi- 
wavelength studies on edge-on discs is equally limited. Two 
objects, HH 30 and the "Butterfly-star" , share a large num- 
ber of features with the disc in CB 26. All three are seen 
circum stellar d iscs, seen almost edge-on. 

Wol f et al] (2003) has compiled a similar data set 
for a model for IRAS 04302+2247, the so called 
"Butterfly-star". The data set includes millimetre maps 
and high-resolution near-infrared images obtained with 
HST/NICMOS. Utilising the same techniques as in this 
study, the authors' main result is that the dust properties 
must be different in the circumstellar disc and in the enve- 
lope. Whilst a grain size distribution with grain radii up to 
100 jam is required to reproduce the millimetre observations 
of the disc, the envelope is dominated by smaller grains 
similar to those of the interstellar medium. This is quite 
in contrast to our model, where we find grains with almost 
ISM-like properties are needed in both, disc and envelope. 
However, the millimetre maps of the model do not suggest 
any presence of a large inner void as the spatial resolution is 
not good enough. In a later investigation (Wolf et al . 2008) 
on the "Butterfly-star" with the Sub- Millimetre Array the 
authors as well discover a "dip" in the horizontal bright- 
ness distribution. According to the model however, in this 
case the effect is to be understood as an effect of the optical 
depth and not as an inner void. 

HH 30 was identified as a circu mstellar disc with a larg e 
inner void present by HH 30 by iGuillote au et al.l (|2008l ) . 
The authors found a inner hole of radius rin = 40±5 AU and 
an outer disc radius of Tout ~ 

128 lb 3 AU. At this numbers 
HH 30 and CB 26 are quite comparable. The separation of 
two central stars is postulated to be about 15 AU, so this 
might also be good guess for CB 26. Yet, in contrast to our 
model for CB 26, no envelope structure is needed to explain 
the appearance of HH 30 in the near-infrared bands. This 
suggests that CB 26 is less dynamically evolved than HH 30. 



Both systems drive a molecular ou tflow, but CB 26 show s 
clear signatures of outflow rotation ("La unhardt et al.]|2008l ), 
while HH 30 does not (Pety et al. 2006j). The reason for this 
difference remains unclear, though the early dynamically 
state of CB 26 as compared to HH 30 may provide the key. 
Yet another difference between HH 30 and our object is 
that the presen ce of cm s ize grains is needed to model the 
value for (3rnm (IWood et a l. 2002). 

As a summary, we see that, although all three ob- 
jects, CB 26, HH 30, and IRAS 04302+2247 feature the 
same structure, a detailed investigation as presented in this 
work shows, that they are systems at different evolutionary 
states. 

6.7. Errors &. caveats 

We conclude this section by discussing simplifications and 
resulting possible sources of errors within the framework of 
our model which have not been mentioned so far. 

A matter that has not been touched so far is the unique- 
ness of or model. Despite the simplifications we applied the 
volume of the parameter space is still to vast to be com- 
pletely scanned in a single study. We therefore employed the 
above described step-by-step technique for the parameters 
to find our best-fit model. 

For instance, we started to model the millimetre maps 
as for those the envelope structure is not a dominant con- 
tribution. Thus, we obtain information about the disc total 
mass as well as its inner and outer radius and the apparent 
absence of large dust grains. The exploration of the near 
infrared images then showed the requirement for a rota- 
tional envelope structure rather then a simple power law 
distribution. By treating the parameters in this indepen- 
dent and serial manner, the order in which they are fitted 
might become an important issue. For instance, the study 
could first have been focusing on the near infrared images 
and wideness of the dust lane. As we know from our expe- 
rience from modelling other objects this wideness usually 
requires small grains at least in the envelope. Modelling the 
SED in the millimetre regime as a second step would still 
have lead us to the overall usage of ISM grains. 

This raises the question, if our model is really a global 
minimum in the parameter space or just a local one and 
if there exists in that space a model, that reproduces the 
observations on CB 26 even better. We exclude this pos- 
sibility. The model we obtained exhibits some unexpected 
features. Therefore, in the framework of our model set we 
thoroughly explored the range of those parameters upon 
which these features sensitively depend. For example, we 
did all the combinations for values of the inner radius and 
the scale height and disc mass. If adjust accordingly, all 
three can produce the plateau observed in the 1.3 mm im- 
age. However, only the larger inner radius prevailed. A 
small scale height might squeeze the dust tight enough for 
high optical depth, but clearly contradicts the wavelength 
dependence of the dark lane in the NIR images. The same 
holds for the maximum grain size. Whereas here we do not 
even have the possibility to mimic ISM behaviour by means 
of disc geometry. 

As a remaining issue, we need to think about parameters 
not varied at all in our modelling effort. As explained in 
sections above, model assumptions such as spherical grains 
versus fractal grains can hold the key to the mass problem 
of the disc. However, this would not alter the model we have 
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at hand, especially this does not provide a hint to the "real" 
global x^-niinimum if one thinks the model is trapped in a 
local minimum. 

This might also be an issue for the stellar parame- 
ters. Varying the luminosity and effective temperature of 
the embedded T Tauri star in general changes the total 
energy throughput in the radiative transfer and the loca- 
tion of the peak of the stellar spectrum in the wavelength 
space. Deviations from our assumed "typical" T Tauri star 
of course will affect the numbers of our best-fit model. For 
instance, since the total mass critically depends on the flux 
in the millimetre regime and this fiux in turn on the total 
energy provided by the central star. Also, the screen in- 
troduced to mimic interstellar extinction is affected by the 
choice of surface temperature. However, no choice of stel- 
lar parameters is able to affect the main conclusions of our 
model. These are the presence of ISM grains in the disc and 
the inner void. 

Despite the discussed caveats, the fact we actually found 
a good model means that we do not need to call for complex 
physics, such as grain growth or dust settling. The data do 
not require this. 

7. Summary 

For a large span of wavelengths, we have compiled a high 
quality data set for the circumstellar disc in the Bok globule 
CB 26. We obtained images in the near infrared and in the 
millimetre regime as well as photometric data and spectra. 
Together with literature values, we have constructed a de- 
tailed model that allows interpretation of observations with 
one single set of parameters. The conclusions we obtained 
are multifarious: 

1 . In order to account for the brightness distribution in the 
1.3 mm map we needed to include an inner hole with a 
radius of 45 AU. Observations to come, especially in 
the sub-millimetre regime are suitable to confirm our 
interpretation of a low optical depth along the line of 
site at 1.3 mm. 

2. Our model very nicely reproduces the prominent chro- 
maticity of the dark lane as seen in the near infrared 
images. 

3. Based on the chosen dust composition (astronomical sil- 
icate, graphite) and the resulting opacity structure of 
the best-fit disc model, we find that the millimetre SED 
indicates that the grains feature the same grain size dis- 
tribution and almost the same upper limit to the grain 
size as the interstellar medium. 

4. The disc is massive with a total dust and gas mass of 
0.3 M© under the assumption of spherical grains and 
Pgrain = 2.5gcm~^, Compared to a mass of the central 
T Tauri star with 0.5 and therefore possibly, but not 
inevitably, unstable. We discussed that due to possible 
fractal structure and other effects the real disc mass may 
be small enough to have a stable system. 
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